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ABSTRACT 

The boundary layer for a swirling flow of incompressible 
fluid in a convergent cone has been investigated. Theoretical 
calculation with boundary conditions more appropriate to 
physically existent situations discount the existence of 
"super velocities" within the boundary layer. Parallel 
experimental investigations demonstrate an interdependence of 
core and boundary layer flows which precludes the maintenance 


Gmeeanitial “ideal flow" conditions. 
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L. INTRODUCTION 

The determination of velocity profiles in the laminar 
boundary layer of a convergent nozzle has been approached 
analytically in several different manners. Taylor [lJ con- 
Sidered the case of a predominant tangential flow with a 
Secondary axial flow. In his analysis, Taylor assumed that 
the tangential flow exceeded the axial flow to the extent 
that the axial component in the inviscid core could be 
neglected. After applying the Pohlhausen method to the 
mementum-integral equations and obtaining solutions for various 
axial positions, Taylor concluded that a condition could 
arise when most of the outflow from the nozzle was fed by a 
boundary-layer current close to the wall of the swirl chamber. 
The radial pressure gradient, which held the particles in their 
Circular path, acted upon the retarded boundary layer to drive 
the fluid along the surface to the exit. 

Wilks [2] developed a similar approach that led to the 
possibility of a 'super-velocity' existing in the boundary 
layer of a convergent conical nozzle. This ‘super-velocity' 
was a velocity within the boundary layer that exceeded the 
free-stream velocity. Wilks considered the flow in the cone 
to be that of an inviscid core with a thin boundary layer. 

The inviscid core had two components of flow - a uniform 
axial flow and a free vortex with its center along the center 
line of the nozzle. 

In this paper an analogous approach was utilized. How- 


ever, the effect of the boundary-layer growth upon the tangential 








and axial velocities in the free stream flow was considered. 
Subsequently, no sSuper-volocities were found to exist. Further- 
more, experimental results prove that full scale effects 


predominate in this problem rather than just boundary-layer 


effects. 
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IIL. ANALYSES 


A. WILKS THEORY 

Using the Wilks approach, a curvilinear coordinate system 
was defined with the x-axis along the wall and the z-axis 
perpendicular to the wall of the nozzle (Figure in Appendix 
A). The shape of the nozzle was determined by the surface of 
revolution of an arbitrary curve r(x). Wilks developed an 
approximate method of solution using the momentum equations 


(Appendix A) 


a 2 | we 2 2 
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aX OZ r-ax dx Yr dx OZ 
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and the corresponding integrated forms of the momentum 


equations 
&  (x{ vu-v)az) + cal (u-U)dz + re (v2-v2) dz (4) 
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where the boundary conditions were 
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The distances in the above equations were non-dimensionalized 
by defining n = z / 6 and X = x / c, where 6 was the local 
boundary-layer thickness and c was the distance along the 
generator to the nozzle to its apex. 

To obtain a solution of the momentum-integral equations, 
it was necessary to assume a velocity profile within the 
boundary layer. A profile that would also satisfy the energy- 
integral equations should produce some improvement over a 
solution to the momentum-integral equation alone. Wieghardt 
[3] developed such a method with the use of two free parameters 
and an eleventh-degree Seneca W. Applying the free para- 
meters and the polynomial to the u-velocity profile and assum- 
ing that the v-velocity profile did not change significantly 
with X, the following profiles were obtained: 

u/U = £1, (n) + Ly*fog(n) + Lo*f3(n) (8) 


where 


f£,(n) = 1- (1 - ny Ol eonne aeen- Gaon) 
fo(n) = (1 - n)8n(1 + 8n + 36n?) 


-(l - n) ®n2(1 + 8n) 


£.,(n) 
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and 

V/V = G(n)e- ez 2n> + n? | (9) 
An evaluation of the axial velocity profile at the wall leads 
to physical connotations for the shape factors, L, and Lo; Clegg 
L, includes the wall shear stress and L» involves the axial 


pressure gradient. 


L,. = ie L. = eo ae 
1 WU 2 2e wv ax 


Using the coefficients Al, A2, A3 and A4 (ratios of character- 
istic thicknesses) as defined in APPENDIX A and introducing 
a term 
2 
UC ic 
Serials, —) (a) = ae) 


the momentum-integral equations became 


r' vu $5,  A2 ey 
eS ce lj) A3 (10) 
x 
q 
Tee U' 6 A4 
Se ee aay 
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As a method of simplifying the solution, Wilks assumed that a 
basically linear relationship existed between the momentum 
thickness (65,) and the boundary-layer thickness (6) and 
between the mixed momentum thickness (Soxy) and the boundary- 
layer thickness. The above equations were then combined 

with the definition of G and its derivative. The following 


relationships were obtained 
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In the preceeding analysis, no restriction was made upon 
the shape of the convergent nozzle. The flow conditions in 
a cone, as well as many other shaped nozzles, do not meet the 
requirements for a Similarity solution of the related boundary 
layer equations. However, in these cases where the profiles 
vary with axial position in a non-linear manner, this present 
approximate method of solution can be Ce oree alee: used. Wilks 
showed that a solution of the similarity equation could be 
developed for any one point along the surface of a cone, but 
that the coefficients of the Falkner-Skan equation would be 
different for each point. For each point, he noted that a 
close agreement existed between the profiles of the approximate 
method of solution and the profiles of the similarity solution. 

For the particular case of a cone, a linear relationship 
existed between the radius and X. The axial velocity as a 
functions of X was then formulated from continuity, and the 
variations in the tangential velocity were derived from con- 
servation of angular momentum. (The details of these functions 


are given in APPENDIX A.) 
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Bs NUMERICAL SOLUTION OF WILKS THEORY 

A solution of the momentum-integral equations would lead 
to the values of the shape factors and thus to the axial 
velocity profile. A step-by-step procedure was used to solve 
the governing equations at various values of X. Wilks 
considered the simplified case where the initial boundary- 
layer thickness on the cone was zero. With this assumption, 
G(0) became zero, which led to L, (0) = 0. Equation 13 then 
meeame a Quadratic im L,(0)- The solution for other values 
of X continued in the following iterative manner: 

(1) Using the values of L, (0) and L, (0), the 
thickness coefficients were calculated. 

(2) The velocities were then evaluated. 

(3) Using equation 13, L, was calculated. 

(4) Applying equation 12, the value of G' was 
determined. 

(5) Using the first two terms of a Taylor series 
expansion, G(X+dX) = G(X) + G'(X) -dXx. 

(6) L5 was calculated from equation 14. 

(7) The values of the thickness coefficients were 
recalculated with the new values of Lj and Ly 
and the procedure from step 2 onward was 
repeated for new values of X. 

In a physical system, a boundary layer could easily exist 
before the flow entered the conical region. The momentum- 
integral equations would then need to be solved for a non-zero 


initial boundary-layer thickness. This was accomplished by a 


i 








Slight modification of the Wilks solution. Using the solution 
ee the quadratic in Tae) fen zereominitial boundary-layer 
thickness as a predicted value of L,(0), an iterative method 
would correct this value for the presence of an initial 
boundary-layer erence eee The method pursued was the following: 
(1) An initial value for L,(0) was assumed, 
(2) The velocities were calculated. 
(3) G and thus L, were determined from the value 
of the initial boundary-layer thickness. 
(4) Using the current values for L, (0) and L, (0), 
the thickness coefficients were calculated. 
(5) The value of L, (0) was determined according 
to equation 13. 
(6) The calculated value of L, was compared with 
the assi-ied value of L,. If the two values 
were not sufficiently close, a value equal to 
the average of the two was assumed and the 
procedure from step 4 onward was repeated. 
If the result was sufficiently accurate, Ly 
and L, were ecaleulated for other values of X 
using the same methods as for a zero initial 
boundary-layer thickness. 
Using the above numerical methods, the momentum-integral 
equations were solved for values of K (swirl factor) of one 
and two. The resultant profiles were compared with the Wilks 


profiles to insure the correctness of the numerical method 
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employed. The solution was then extended for values of K up 
to six with various initial thicknesses and initial velocities. 
Several profiles that were obtained in this manner are shown 
in Figures 1-7. The solutions indicated that the degree of 
super-velocity increased dramatically with K and X. The 
value of the momentum thickness became negative under many 
circumstances, and for high degrees of super-velocity, the 
displacement thickness also became negative. Under all of 
these conditions the boundary-layer growth continued to 
increase with X. The numerical analysis deteriorated for 
values of X much greater than 0.3 in most of these solutions. 
There was also a failure in the numerical solution for ratios 
Seetangential velocity to axial velocity in the range of three. 
It was believed that these problems were caused by the numer- 
ical method employed and were not inherent problems of the 
analytical model. The quantity (A2 - A4) could be very close 
to zero. Since one of the terms in the calculation of Ly 
required division by ths above quantity, overflow in the 
computer could easily be caused when this quantity approached 
zero. A different final arrangement of the momentum-integral 
equations would avoid this problem and allow for a solution 


at higher values of X. 
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FIGURE 1. VARIATIONS IN PROFILE WITH K. xX=0.1 G(0)=0 


(a) K=0 , (b) K=l1 , (c) K=2 , (d) K=3 , (e) K=4 
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FIGURE 2. VARIATIONS IN PROFILE WITH K. X=0.1 G(0)=0.52 


(a) K=l , (b) K=2 , (c) K=3 , (da) K=4 
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FIGURE 3. VARIATIONS IN PROFILE WITH K. X=0.2 G(0)=0.52 
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FIGURE 4. VARIATIONS IN PROFILE WITH X. K=1 G(0)=0.52 
(a) X=0.00 , (b) X=0.05 , (c) X=0.10, (dad) X=0.15 


(e) X=0.20 , (£) xX=0.25 
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FIGURE 5. VARIATIONS IN PROFILE WITH X. K=2 G(0)=0.52 
(a) X=0.00 , (b) X=0.05 , (c) X=0.10 , (d) X=0.15 


(e) X=0.20 , (£) X=0.25 


Za 








FIGURE 6. VARIATIONS IN PROFILE WITH X. K=4 G(0)=0.52 
(a) %=,00 , (6) xX=.05 , (c) X=.10 , (dad) X=.15 , (e) tee 
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FIGURE 7. 





VARIATIONS IN PROFILE WITH G. K=2 X=0.2 


(a) 6(0)=0.00 , (b) 6(0)=0.10 , (c) 6(0)=0.15 
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om PRESENT THEORY 

The Wilks approach to the momentum equation assumed that 
the tangential velocity at the outer edge of the boundary 
layer could be described as A/r. This was only appropriate 
if the boundary-layer thickness was very small compared to 
the value of r. The numerical solution of the Wilks method 
led to a boundary-layer growth that could be of the same 
order as r for flow in air where r, was about one inch and the 
initial axial velocity was one foot per second. In these 
cases, the radius minus the boundary-layer thickness should 
be used to determine the tangential velocity at the outer 
edge of the axial boundary layer rather than using the radius 
alone. Thus V = A/(r - 46). 

The introduction of this new value of V into the boundary 
conditions of the momentum equations led to a different form 
of the v momentum-integral equation. The derivation of the 
equation is shown in APPENDIX B. This approach required the 
definition of a new thickness coefficient A5 (APPENDIX B). 

The v momentum-integral equation then became 

Pepe hese +n (F- Fe 
In the above equation, Wilks assumption of a basically linear 
relationship between the characteristic thicknesses was used. 
miemwcerms in the braces are added to the Wilks v momentum-— 
integral equation by the new Seer cice vix., V = A/(r - 4). 

If the boundary layer is thick enough to effect the 


boundary condition for the tangential flow, the calculation 
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of the axial velocity should also take this boundary-layer 
fmewth into cons#deratzon. In using@the conditrons of conrer— 
nuity of flow, the radius minus the displacement thickness 
should be used rather than the radius alone. The displacement 
muuckness, by its definition, is the distance the potential 
flow has been displaced by the presence of the boundary layer. 
The effect of this consideration would be to increase the 


axial core velocity as the displacement thickness grew. 


1. NUMERICAL SOLUTION OF PRESENT THEORY 

In order to calculate U and its derivative, U', the value 
of the displacement thickness and its derivative were needed. 
The ratio of the displacement thickness to the boundary-layer 
thickness was one of the integrals used in the calculation 
of the thickness coefficients. The boundary-layer thickness 
at any point could easily be calculated from the value of 
G once the kinematic viscosity and length of the nozzle to 
its apex had been specified. The calculation of the derivative 
of the displacement thickness was not as direct. Wilks based 
his solution on the assumption that a baSically linear relation- 
ship existed between the momentum thickness and the boundary- 
layer thickness. An extension of this assumption would be 
to presume a baSically linear relationship also existed between 
the displacement thickness and the boundary-layer thickness. 
Thus the ratio of the displacement thickness to itS derivative 
should be equal to the ratio of the boundary-layer thickness 


mo LtS derivative. 
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With the new equation for the v momentum integral, the 


numerical solution proceeded in the following manner: 


(1) 


(2) 


(3) 


(4) 


(5) 


Using the initial boundary-layer thickness, 

the values of G(0) and i510) were icaleulacces 

A predicted value for L, (0) was assumed from 
the Wilks solution for zero initial boundary 
layer thickness. 

The distance along a flat plate that would be 
reguired to produce a boundary-layer thickness 
of the same size as the initial boundary-layer 
thickness was calculated. Using this distance, 
the derivative of the boundary-layer thickness 
on a flat plate at this point was calculated. 
From this estimate of the derivative and 
assuming an initial value for Dl from the Wilks 
solution, the initial value of the derivative 
of U was evaluated. 

Usage ene cueneht yates On ye clCm omen c 
values of the thickness coefficients were 
calculated. 


The value of 6'/dé was calculated, where 


or A4 Us Aa AS V 
oe SG Oi tc nem maven 
os as . 
as) 5 a (16) 
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(6) 


(7) 


(8) 


Co 


(10) 


The new value of Ly was caleulated cdizom 


; ? 2 
Ly = p2-c- | O (2 ape eee re a 
U xv U2 é 
7) 


The calculated value of L, (0) was compared with 
the assumed value. If the two values were not 
sufficiently close, the thickness coefficients 
and velocities were re-evaluated until the 
values of L, (0) for two consecutive iterations 
were sufficiently close. 


The value of G' was determined from 

G' =G a + 

The values of G and Ly for the next point were 
then calculated. 

The step-by-step evaluation of the momentum- 
integral equations was continued for new values 


of X by returning to step 4 and remembering 


that step 7 1s used only when X=0. 


U and V as functions of X are given in APPENDIX B. 

Using this numerical technique, the momentum-integral 
equations were solved for various values of initial boundary- 
layer thicknesses, initial eaieeecey and K. Since this 
method required the division by G in the calculation of equa- 
tion 16, the analysis was not valid for the case where the 


initial boundary-layer thickness was zero. However, for very 
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small values of G other than zero, there was no restriction 

on the initial boundary-layer thickness. The solution of 

these equations showed no sign of super-velocity. In several 
cases the step-by-step evaluation continued smoothly for values 
Sex up to 0.6. The boundary-layer growth reached a maximum 
value when X was in the range of 0.25 and then decreased with 
X. This numerical analysis appeared to lose its validity when 
the momentum thickness became negative, which occurred at a 
much higher value of X then in the Wilks solution. At this 
value of X where the momentum thickness became negative, dis- 
continuities appeared in the boundary-layer growth. Figures 8-12 
show a comparison between the profiles of the Wilks solution, 
aS modified for initial boundary-layer thickness, and the 


profiles obtained by the above method. 


aie DISCUSSION OF THEORETICAL RESULTS 

Wilks based his solution of the flow in the boundary layer 
Of a convergent nozzle on several assumptions. He first con- 
Sidered the inviscid core as a free vortex with uniform axial 
flox. It was then assumed that the boundary layer would be 
thin enough to prevent any effect of its growth upon the core 
velocities. Thus the core velocities were dependent only upon 
the geometric configuration of the convergent nozzle. Wilks 
also assumed that a linear relationship existed between the 
various characteristic thicknesses of the boundary layer. With 
these assumptions, the solution of the momentum-integral 
equations led to the existence of a super-velocity in the 


boundary layer. 
ao 





For low velocity air flow (required to maintain the 
laminar conditions) the boundary-layer growth was not negli- 
gible. By including this growth in the calculations of the 
boundary conditions, the solution no longer led to super- 
velocities. Increased swirl only Beene the ratio of the 
displacement thickness to the boundary-layer thickness 


eeeughtily. 
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FIGURE 8. COMPARISON OF THEORETICAL CURVES. K=l x=.1 
G(0)=0.52 6 (0) /r,=0.130 (a) Wilks theory , (b) present 


theory 
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FIGURE 9. COMPARISON OF THEORETICAL CURVES. K=1l X=0.2 G=0.52 


6 (0) /r=0.130 (a) Wilks theory , (b) present theory 
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FIGURE 10. COMPARISON OF THEORETICAL CURVES. K=2 X=0.1 
G(0)=0.52 S0)/E5=0s130 (a) Wilks Cheer, |; (b)epresent 


theory 
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FIGURE 11. COMPARISON OF THEORETICAL CURVES. K=2 X=0.2 
G(0)=0.52 6 (0) /x=0.130 (a) Wilks theory , (b) present 


theory 
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FIGURE 12. COMPARISON OF THEORETICAL CURVES. K=4 X=0.1 
G(0)=0.26 G7 ae Uae (a) Wilks theory , (b) present 


theory 
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fil: EXPERIMENTATION 


A. DESIGN OF APPARATUS 

The analytical approach assumed that the flow entering 
the convergent nozzle was a combination of a free line-vortex 
and uniform axial flow. To test the theory, an apparatus was 
required to approach this assumed velocity distribution. There 
are several methods of introducing a fluid into a chamber and 
causing a Swirling motion. However, this tangential velocity 
would have to approximate the conditions of a free corte that 
is the tangential velocity should be inversely proportional to 
the .Yadius. Several methods were considered, including rotating 
a section of the swirl-producing chamber, drawing the fluid into 
the chamber from a stagnation condition while using vanes to 
create the tangential motion, and rotating blades inside the 
chamber. The approach chosen was to introduce the fluid through 
one tangential port with variable cross-section. In tests of a 
Hilsch tube, Lay [4] used a single circular port to introduce 
a jet into the chamber. This was the only source of fluid in 
his tests, and he reported the velocity profiles that were 
produced in this manner. However, he did not indicate the 
radial profile of the velocities in a region close to the wall. 
The results indicated that the tangential flow might closely 
approximate the conditions for a free vortex over most of the 
radius, departing from these conditions in the center portion. 
The distribution of the axial velocity in Lay's tests suggested 
that the superposition of pipe flow upon this flow could result 


in a nearly uniform axial flow. 
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The source of fluid for this apparatus was a compressed 
air tank. To insure laminar pipe flow of the air at speeds 
of at least two feet »er second, a pipe diameter of two inches 
was chosen. Axial pipe-flow could be produced prior to entry 
of the axial flow into the swirl chamber. The shape of this 
axial flow could be varied by adjusting the axial position of 
a honeycomb flow-straightener. If profiles different from 
pipe flow were required to achieve the uniform axial profile 
at the entrance of the cone, the flow could be adjusted by 
only pasSing air through selected annular portions of the 
honeycomb flow-straightener. 

The tangential flow was introduced through a rectangular 
slot in the wall of the swirl generator. One wall of this 
Slot was tangent to the inner surface of the chamber. The slot 
was 0.25 inches high and the width could be varied from 0.3 
to 1.25 inches. This variable cross-section permitted an 
adjustment in the ratio of the velocity to the quantity of mass 
introduced by the tangential port. This flexibility in the 
system was added to insure that the proper profiles and values 
of K could be obtained. Tests were then performed to determine 
the axial distance required for the development of the appro- 
priate flow conditions. These preliminary tests indicated that 
a range between one and two feet was required for the development 
of the tangential jet into a free vortex. 

Lay's tests indicated that there would be a region close 


to the wall where the tangential flow wouid be retarded by 
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viscosity. An axial combination of pipe flow and the tangen- 
tially introduced flow would also have a region close to the 
wall where the conditions for uniform axial flow would not be 
met due to the viscous forces present. With these two boundary- 
layer growths under consideration, ine outer annular portion 
of the flow was bled off, with only the "flat" inner portion 
of the axial flow entering the cone. This bleed-off was 
accomplished by placing a cone of initial radius 0.77 inches 
at the end of the tube which had an inside radius of one 

inch. The position of the cone relative to the end of the 
tube could be varied in order to alter the amount of flow bled 
off into the atmosphere. The entry region into the cone 
consisted of a one-inch straight section. The purpose of this 
section was to insure the continuation of the proper flow 
characteristics after bleed-off. (See Figure 25) 

The cone, a piece of machined plexiglass, featured a 
semi-angle of 4.4 degrees. The distance along the generator 
of the cone to the apex, termed c in the analysis, was ten 
inches. The truncated cone used in the test apparatus extended 
over 0.35 c. Thus, tests could easily be performed over the 
region where X varied from 0.0 to 0.3. Two hot-wire probes 
were introduced from the exhaust end of the truncated cone. 
One configuration was sensitive to axial and radial components 
alone. The other probe could detect tangential and radial 
velocities. It was assumed that the radial component of flow 


was negligible, especially close to the wall. The traversing 
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mechanism moved in two directions, allowing for motions in 


both the x-direction and the z direction. 


Bs EXPERIMENTAL RESULTS 

The development of flow conditions that would approximate 
the desired conditions was basically a trial-and-error pro- 
cedure. By adjusting the slot width, the development length 
for the tangential flow (by changing the length of pipe between 
the swirl generator and the nozzle), the bleed-off level, and 
the radial distribution of the axially introduced flow, it was 
Pe-ssitble to achieve several flow conditions that approximated 
the analytical boundary conditions at the entrance of the cone. 
However, in all the cases tested, the velocity distribution in 
the inviscid core chan ,ed drastically with X. The radial dis- 
tribution of the axial and tangential velocities for several 
values of X with swirl ratios (K) of 0.6, 0.8, and 0.95 are 
Shown in Figures 13-24. These profiles indicate several 
different flow relationships. For the case where K=0.6, the 
degree of swirl had very little effect upon the axial flow as 
it moved through the cone. The uniform axial core was gradually 
replaced by a flow that resembled the distribution in pipe flow, 
with most of the flow occuring in the central portion, and the 
perimeter retarded by viscous forces. With the slight increase 
to K=0.8, there was a great change in the radial distribution 
of the axial flow. In this case, the free vortex seemed to be 
maintained for all test stations. At X=0.3, this tangential 
motion had thrown a great deal of the fluid into the peripheral 
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portion of the cone, leading to a High axial yeltocimancam ene 
wall. oveaee this high velocity was still in the core 

(z=0.2 inches) and not in the boundary layer. In the tests 
where K=0.95, the distance over which the tangential flow 
eepeared as a free vortex was less then in the previous cases. 
This smaller distance resulted in several changes in the 
tangential velocity with X. As a result of this changing 
tangential profile, there were several in the axial flow. In 
sections where the tangential flow appeared as a forced vortex, 
there was one form of axial flow, where high speed occurred 
away from the wall, while in annular sections where the tan- 
gential flow resembled a free vortex, the high speed axial flow 


occurred nearer the wall. 


Cc. DISCUSSION OF EXPERIMENTAL RESULTS 

In both of the analytical approaches, the inviscid core 
was considered a combination of uniform axial flow and a free 
vortex. In the test cases, the axial flow at the entrance of 
the nozzle was within five per cent of uniform axial flow. The 
tangential flow approximated a free vortex over most of the 
radius in the entrance region. However, this flow could not 
be maintained down the cone. Swirl induced a high velocity 
axial flow near the wall when the tangential flow was in the 
form of a free vortex. This region of high axial-velocity flow 
was far enough away from the wall (a distance greater than 0.15 


inches) to call it part of the inviscid core. Thus swirl does 
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andeed cause changes in the axial flow, but these jenangessecems 


in the core and not in the boundary layer as Wilks suggested. 
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FIGURE 13. EXPERIMENTAL PROFILES K=0.6 X=0.0 
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FIGURE 14. EXPERIMENTAL PROFILES. K=0.6 xX=0.1 
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FIGURE 15. EXPERIMENTAL RESULTS. K=0.6 X=0.2 
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FIGURE 16. EXPERIMENTAL PROFILES. K=0.6 xX=0.3 
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FIGURE 17. EXPERIMENTAL PROFILES. K=0.8 xX=0.0 
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FIGURE 18. EXPERIMENTAL PROFILES. K=0.8 xX=0.1 
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FIGURE 19. EXPERIMENTAL PROFILES K=0.8 xX=0.2 
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FIGURE 20. EXPERIMENTAL PROFILES. kK=0.8 X=0.3 
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FIGURE 21. 


EXPERIMENTAL PROFILES. 
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K=0.95 





X=0.0 
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FIGURE 22. EXPERIMENTAL PROFILES. K=0.95 xX=0.1 
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FIGURE 23. EXPERIMENTAL PROFILES. K=0.95 X+0.2 
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FIGURE 24. EXPERIMENTAL PROFILES. K=0.95 xX=0.3 


53 





*“NOILWYNDSIANOOD NOILOUS LSAL °*Se dunold 


uoTbsrz JJO-psslq i> 


/ eee eee 


uOoOTRODS 3S9R 
uotbsr yusuedoToeaaep 


uotThsr souerzQusa., 


54 





IV. CONCLUSIONS AND RECOMMENDATIONS 

Previous work in this area considered the effect of 
Swirling flow upon the velocity profiles in the boundary 
layer. In this paper, the effect of the boundary layer growth 
upon the core velocities was utilized in the analysis. How- 
ever, the experimental tests showed that this consideration 
was not enough. The swirl also affected the axial core flow, 
and the resultant changes are much more drastic than any varia- 
tion in the boundary layer flow. Rather than applying the 
Navier-Stokes equations to the boundary layer alone, it would 
be much more appropriate to apply them to the whole flow field. 
In addition, the developed theory for the boundary layer assumed 
that a basically linear relationship existed between the 
various characteristic thicknesses. Since these thicknesses 
are functions of the shape factors, which are in turn functions 
of X, it would seem mandatory to extend this analysis without 
these assumptions of linear relationships. However, in order to 
complete such an analysis, more information needs to be known 


about the effect of swirl upon the core velocities. 
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APPENDIX A 
DERIVATION OF EQUATIONS ACCORDING TO WILKS 
The boundary of the convergent nozzle is formed by sur- 
face of revoltuion of a curve r(x), where x is a distance 
measured along the wall. The distance from the wall is az. 
The potential flow has the two components, V and U as shown 


below 





Acceleration in the orthogonal curvilinear coordinates is 


defined as 


ie 


= Do 9 = 
a= 5H = =2 + WV (1/2q7) -qx® 


In the st direction this becomes 


1 1 1 2 2 3 3 
Z Z 2 
; Jt qy (2) qd. (~ : qd, (~*3) 
Oe hy ou. No ous h3 ous 





af y= Xe eee, U3 


i 
N 


116.4. Wy eee VOR oe 
hy = 1} ’ ho soe I h3 = 


the acceleration can be written 


_ ou... ef v ou CU Os 

oe at || esa SO ae eee: 
ow Ow ov Ow aw ve ax 

44 = ost tu Ox ty 96 +twWQRZ7 YF Dz 
BY ie ew aa ON, pov... Wau os 
=o Ot U Dx r 00 ne OZ t r Ox. 


For axisymmetric flow, all of the partial derivatives with 
respect to the angular displacement are zero. Since velocities 
in the direction perpendicular to the wall (the z direction) 
may be considered small compared to the other velocities, 

the momentum equation in the z direction may be neglected. By 
an order of magnitude analysis, the boundary-layer equations 


for steady laminar flow with constant properties become 


é 2 
ou +w2Riv Slit Py ye au 
aX az jae CL p ox JZ 

2 
Frere oe a oc he 2 oy 
ox. OZ r dx 0 08 OZ 
on (r-u) + 2 (rew) = 0 
OX OZ 
with the boundary conditions 
u=v= 0 on z= 0 
u> U(x) ard” v-> V(x) = A/r(x) as zs © 
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With the above boundary conditions, sine pressure sgradlcnmuser 
the x-direction can then be written 
0 


eS du y2 xt 
ox & ax Y r 


Di 
rg 


Integrating the u-momentum equation across the boundary layer 


gives 
(a. = Pe a | 
dU Clee Ve Cie au os __., fou 
J (uo Bewe r dx ue + v2 E \az =v (Re 


from continuity 


fut 


| Z 
w=-- | Ce (reu) dz 


° 
os 


Thus 


| (» ou ~ i ou o (reu) az - v2 re’ =a au 
0 ie 9 Ox ie 


ox OZ dx 
pen —_ su 
+ v2 E)az = »( 32 
O 
The second term 
heen ae 
, 22] ee <2 
o = OZ a ox 


can be integrated by parts. 
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00 Z. 
m p | ~ | p dm=u | <— (xu) dz ~ | u a (r-sWjiazZ 
X ox 
O O O O O 
= U | 9 (reu) dz - | u 3_ (reu) dz 
0 OX 0 om 


The u momentum-integral equation then becomes 


| (u ou = Ue eu) es a (r°u) - U -) az 
0 OX Yr Ox Y ox dx 


fe | (v2 ~ v2) r’ a7 Sy) ( $2} 
O x O 


The first integral can be written 
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The integral form of the v-momentum equation is 


| uv ov 4 a ov + usv x=} dz-=- ) (33) 
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imeroducing the Gentinus ty sequdeien, cht s becomes 


+e @ a @ 1 


ec 


dv _ 1 ov an - Wien ar oa OV 
[es aa ce O32 qe) a = » (32) 


Performing an integration of the second term by parts 
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The left hand side of the momentum-integral equation may now 


bemwritten as 











. ov _ V 3 ‘ wd ; “uev dr 
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Thus the v momentum-integral equation is 


with 


e. :? | 
r ax O 
the 
v = 0 
97y = QO 
927 
ve = Vion) 
OV 
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The u momentum-integral equation is 


Rie 


a | > 
Mae 


|» 


+ 


en, 


He 
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le elt 

ou 

Aol 0 when 2 > @ 
ie ' 


In order to obtain a solution to the momentum-integral 
equations using the approximate method of von Karman, it is 
necessary to assume suitable velocity profiles. A direct 


analogy with the Pohlhausen method would lead to 


u/U = f£f(n,L}) V/V = g(n,Lo) 
where n = 2Z2/6(x). The continuity equation leads to a stream 
function 


Ww = rvreUrh(x,n) 


The axial velocity can then be written 


_ Oo 
u = sy bh Gx,n) 


If the tangential velocity is of the form v = Veg(x,n) , 
the momentum equations may be written 

h"' + acheh" + be (1-h'?) + d*(1-g2) = 0 

6 eaeh=g* = 0 

Where a, b, and d are constants as defined by Wilks. Boundary 
conditions require that g" evaluate at the wall is zero. A 
differentiated form of the v-momentum equation gives 

an e+ aimee. “aireg ) = ~0 
baumee both h’ and g” evaluated at the wall are zero, it is 
apparent the g"* evaluated at n = 0 must also be zero. This 


would require the second and third derivatives of g to be 
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itndependent of x. With this understanding, the two parameters 
can be used to obtain the axial velocity profile while the 
tangential profile is considered dependent upon n alone. 

Thus 

a and v/V = g(n) 

Introducing a dimensionless length X=x/c (c being the length 


u/U = £(n,L 


of the cone along the generator), the u momentum-integral 


equation becomes 


1 a eal 
aL a 1 dU 
Bae ax | r+U2+6 | £¢é-10¢8] + = ax U°d | Cia alia 
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O O 
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with the boundary conditions 
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The v momentum-integral becomes 
a 2 1 Vv 
pee | eee TO f(l-cq) on = — V-g' (0) 
err ax O 6 


with the boundary conditions 


Ci 
on n= 0 


g 0 
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g' = 0 on n=l 


g" 0 
Here it was assumed that the boundary conditions could be 
transferred from infinity to the boundary-layer thickness, 
that is where 0 without appreciable error. 

In an attempt to obtain an improvement in the basic 
approximate method of oeonp a profile can be obtained that 
will also satisfy the requirements of the energy-integral 
equation. Weighardt [3] developed such a method uSing two- 
parameter profiles. This profile can be represented with the 


aid of an eleventh order polynomial. Thus 


uf/U = £,(n) + Li °f. (n) + L,°f.(n) 
where | 

£,(n) = 1- (1-n)8(1 + 8n + 36n2 + 120n3) 

f(n) = (1-n)®n (1 + 8n + 36n%) 

f,(n) = -(1-n) ®n2(1 + 8n) 


Evaluating the functions at the wall (n=0) yields 


au) . wu oe) uy 
(32), (ae), pea! 
thus 
- ota 
1 uU 
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Also 


VU U sau va dx 
Se oe eae eel at n=0 
6° Cwicix Ce eX 


where £"(0) = ~2 Lo 


This leads to 
| 2 


ioe = 
ZC U U ae 


The integrals in the momentum integral equations can be 
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ee 
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expressed as various characteristic thicknesses. 


1 
64 = 6 | (1-f£) dn displacement thickness in the 
O x direction 
ora 
Oo, = § | f (1-f£) dn momentum thickness in the x 
O direction 
1 
S oxy = 6 | (1-g) dn mixed momentum thickness 
O 


mame the definition of a new thickness 


- 2 
$5, = 8 | (1-g2) an 
O 


the momentum-integral equations can now be written 
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Differentiating the appropriate terms, these equations 


Simplify to 
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Pamencegration of the polynomials in the above integrals led 
to 
Dl = 0.3333333 - 0.0454545 L, + 0.006060 L, 
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D2 0.0740342 + 0. 01491908) — 0.001 eSe3 L 


Hl 2 


II 


2 


- 0.0045481 te ~ 0.00008286 L, 


Tomo L203 LjL5 


D3 


0.08461538 + 0.02435846 L = 0.00 293040 L, 


D4 0.4174603 


ll 


Wilks made the assumption that a basically linear relationship 
existed between the boundary-layer thickness and the momentum 
thickness. This assumes that D2 does not vary significantly 
with X. A basically linear relationship between the mixed- 
momentum thickness and the boundary-layer thickness was also 
assumed. 

With the introduction of a dimensionless parameter 


vs 
V 


Q| am 


G = 


the foregoing equations can now be written 


2 
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In addition to the above relationships, the previously derived 


equation for the shape factor L. is also important, viz., 


2 
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1, = Sie y) ri 
Z cate) U 1G 
The following characteristics must be considered for the 
perfectly conical nozzle. 
= - X 
x (X) rol ) 
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Conservation of angular momentum requires that 
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APPENDIX B 
DERIVATION OF PRESENT THEORY 

The concept 6f a free vortex in the re etal core requires 
that the tangential velocity be inversely proportional to the 
radius. In the case where boundary layer growth is large, the 
boundary layer thickness can no longer be considered thin in 
comparison to the radius. With these considerations, the 
boundary conditions for the momentum equation become 

v > V(x) = A /(r-4) as Z2> ® 
This new definition of V(x) has no effect upon the development 
of the u momentum-integral equation. However, there is an 
effect in the v momentum-integral equation. The development 


is the same as described in Appendix A up to the point 
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OX r=—6 ax ax 


Thus the momentum-integral equation becomes 
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Simplifing 
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Introducing the velocity profiles as in Appendix A, the momentum- 


integral equations and the appropriate boundary conditions 


become 
| | 1 1 
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ene) 6 a =e 
(iv) Sa) on n= 0 
(v) ai = 0 


Again, it was assumed that a basically linear relationship 
existed between the boundary-layer thickness and both the 
momentum and mixed-momentum thicknesses. Defining the 
coefficients Al, A2, A3, and A4 aS in Appendix A, and intro- 


ducing a new term 


BO 2.0 ae ee 


D3 
| £ (oie 


the momentum-integral equations then becam: 
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The second equation can be rewritten 
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Thus, the shape factor Ly can be obtained from 
ma 8 2 ' 
a2 = o[2° (2+ ary + Ea - Yas) + 2) 
U iF U 6 


and the new values for G and Lo can be calculated from the 


value of G', where G' is now 


eal 





and 


As the boundary layer grows, the axial velocity as defined 
in APPENDIX A for a cone will no longer be valid. The increas- 
ing boundary-layer thickness will cause the flow in the axial 
direction to accelerate more than would be caused by the 
convergence of the cone alone. The concept of displacement 
thickness would take into consideration this change in the 


effective radius of the cone. Continuity now yields 


v4" [r, - 81 (0)] 2m : 0x) +[x5: (1-X) - 6,00] om 
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With the new definition of V(X), the conservation of angular 


momentum yields 
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